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The understanding of transport processes inside the stars is one of the main goals of asteroseismology. Chemical turbulent 
mixing can affect the internal distribution of /j, near the energy generating core, having an effect on the evolutionary 
tracks similar to that of overshooting. This mixing leads to a smoother chemical composition profile near the edge of the 
convective core, which is reflected in the behavior of the buoyancy frequency and, therefore, in the frequencies of gravity 
modes. We describe the effects of convective overshooting and turbulent mixing on the frequencies of gravity modes in 
B-type main sequence stars. In particular, the cases of p-g mixed modes in Cep stars and high-order modes in SPBs are 
considered. 
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1 Introduction 

The standard stellar e volution modeling in cludes only the 
mixing-length theory ( Bohm- Vitensdl 1958 ) to describe the 
convective mixing, and, for low mass stars, also the effect 
of microscopic diffusion. Contrarily to predictions of this 
modeling, some chemical transport must also occur in the 
stellar radiative regions if we want to account for some ob- 
servational facts. In particular, the comparison between the- 
oretical models and obs ervational data from binaries and 
stellar clusters (see e.g. lAndersen et al. Il990t iRibas etalj 
2000) has shown that the standard model underestimates 
the mixing in the central regions of the stars. It is largely 
accepted that some extra-mixing takes place just at the bor- 
der of the convective core that develops for stellar models 
with masses larger than ~ 1.2 M Q , but there is not consen- 
sus about the physical process es responsible for th is mixing: 
convective overs hooting (e.g 
scopic diffusion (Michaud et al 



Schall er et al.lll992l) . micro- 
2004), r otationally ind uced 



mixing (see e.g. lMevnet & Maeder 2000t lMaederl2003 , and 
references therein), mixing by non-adiabatic propagation of 
gravity waves excited a t the boundary of the convective core 
dYoung & Ariiettll2005l) . etc. 

Generally, stellar modeling includes an overshooting pa- 
rameter (aov) to account for the extra-mixing above the 
limit of the convective core, whose value has been estimated 
by fitting of stellar clusters or binary systems. Recently, as- 
teroseismic determinations of the overshooting parameter 
have also been done by using the advantageous properties of 
low order p and g modes detected in some Cep targets (see 
e.g.lPamvatnykh et alfcOollAerts et al.l2003llBriquet et al 



19911 lAudard & Provostl[l994l) . The shape of the composi- 
tion transition zone is also a matter of great importance as 
far as asteroseismology is concerned. In particular it signif- 
icantly affects the term V M appearing in the Brunt-Vaisala 
frequency and plays a critical role in the phenomenon of 
mode trapping. What w e show in this co ntribution is that 
(as already suggested by IGoupil & Talonll2002l for S Scuti 



stars) the oscillation spectrum of B-type stars is also very 
sensitive to the slope of the chemical gradient in the cen- 
tral region of the star, and therefore, to the kind of transport 
processes that produced it. 

The models with overshooting and turbulent mixing we 
computed to analyze how different extra-mixing processes 
affect the oscillation frequencies are described in section|2] 
and the effects of these transport processes on the stellar 
structure and on the spectrum are reported in section [3] In 
a first approach we consider a parametric turbulent mixing 
without caring about its physical origin. If this chemical 
transport were generated by rotation accordi ng to the for- 
malis m implemented in Geneva code (see. e.g. lEggenberger et aL[ 



2007), for instance, the effects due to the angular momen- 



tum transport should also be taken into account in the os- 
cillation spectrum. In section 13.31 we present an estimation 
of these effects for a typical (3 Cep star and in section 0] we 
give our conclusions. 

2 Stellar models 

To study the effects of different extra-mixing processes in 
SPB and (3 Cep pulsators we focus on models of 6 and 
10 M Q respectively. As a first approach, the turbulent mix- 
ing has been modeled by a diffusion process with a paramet- 
crossing dAizenman et al J 1 977b are very sensitive to the ex- ric turbulent diffusion coefficient D T that is constant inside 
tension of the mixed region (e.g. lDziembowski & Pamvatnvk hf he star and independent of age. This is a very rough approx- 

imation without any a priori physical justification. We com- 

* e-mail: j.montalban@ulg.ac.be puted also some models using the Geneva code that include 



2007al) . I n fact, the frequency o f modes close to an avoided 
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Fig. 1 HR diagram showing main sequence evolution- 
ary tracks of 10 M Q models computed with CLES evolution 
code with an overshooting parameter «ov = 0.1 (dashed 
line), and with a constant turbulent diffusion coefficient 
Dt = 7 x 10 4 cm s~ 2 (solid line). The dotted line cor- 
responds to the evolutionary track of a 10 M Q model com- 
puted with the Geneva evolution code assuming an initial 
rotational velocity V IQ t — 50 km s — 1 . The symbols in the 
middle of MS track indicate the location of models with 
X c = 0.3. 



the treatment of rotation and linked transport processes such 
as described in Maeder & Zahn (1998) and Maeder (2003). 
The comparison between both series of models shows that 
chemical composition profiles in the central regions pro- 
vided by a uniform and time independent diffusion coef- 
ficient represent a first approximation, at least for massive 
stars, of the effect would be produced by such a rotationally 
induced chemical transport. 

The models with overshooting wer e computed with C LES1 
(Code Liegoise d'Evolution Stellaire IScuflaire et alJ l2007). 
In this code the thickness of the overshooting layer Aov 
is parameterized in terms of the local pressure scale height 
H p : A v = «ov x (min(r cc , H p (r cc )), where r cc is the 
radius of the convective core given by the Schwarzschild 
criterion and aov is a f ree parameter. The values of aov 
used in this study were 0, 0. 1 and 0.2. 

The values of the parameter Dt were chosen in order 
to be close to the value of the chemical diffusion coefficient 
near the core provided by the Geneva models. SPBs and 
j3 Cep pulsators are considered as slow or moderate rota- 
tors. SPBs show a ty pical rotational velocity of 25 km s _1 
( Briquet et al. 2007bh . whereas the range of projected rota- 
tional velocity in (3 Cep stars ext ends from to 300 km s" 1 
with an average of 100 km s" 1 dStankov & Handlerll2005l) . 
The Geneva code calculations for 10 M Q models provide 



values of the chemical diffusion coefficient near the con- 
vective core with Xq = 0.3, of the order of 5 x 10 4 cm 2 s _1 
for an initial rotational velocity (Vi) of 20 km s -1 , 7 x 
10 4 cm 2 s _1 for 14=50 km s -1 , and 1.6 x 10 5 cm 2 s _1 for 
Vi=100 km s _1 . On the other hand, the effect of an initial 
rotational velocity of 25km s _1 on the central hydrogen dis- 
tribution of a 6 M© model is well mimicked by a Dt ~ 
5000 cm 2 s _l . The results presented in this paper concern 
mainly the parametric models with Dt = 5 x 10 3 cm 2 s 
for SPB model and D T = 7 x 10 4 cm 2 s" 1 for (3 Cep 
one, and they will be compared with overshooting models 
closely located in the HR diagram, that means aov=0 for 
the SPB case and 0. 1 for the (3 Cep one. 



3 Overshooting versus turbulent diffusion 

3.1 Effects on the structure 

Whatever the origin of the chemical extra-mixing at the bor- 
der of the convective core, the effect on the evolutionary 
track is an increase of luminosity and of the duration of the 
main-sequence phase. The evolutionary tracks of 10 M© 
models with different treatments of chemical mixing are 
shown in Fig. [2] Though the HRD location of turbulent mix- 
ing models is quite close to that of models computed with 
overshooting, the chemical composition gradient near the 
core is significantly different (see Fig. [2] upper panel) and 
therefore also the properties of the Brunt-Vaisala frequency 
(Fig. [2]lower panel). Note that the hydrogen abundance pro- 
file near the core that results from Geneva code (including 
a consistent treatment of rotationally induced transport pro- 
cesses) shows a behavior similar to that obtained by using 
the parameterized turbulent mixing in CLES. The values 
of the chemical diffusion coefficient induced by rotation in 
Geneva models for three different evolutionary stages are 
shown in Fig. 13. II for 10 M© models with an initial rota- 
tional velocity Vi = 50 km s _1 . We note the slight tempo- 
ral variation of this diffusion coefficient near the convective 
core boundary (ra/M ~ 0.15 — 0.25) and the fact that the 
value of Dt we chose (horizontal line) rather corresponds 
to an average of Geneva model values. 



3.2 Effect on oscillation spectrum 

SPB and (3 Cep stars show oscillation modes excited by the 

classical ^-me chanism in the iron group opacity bump at 

T ~200000 K dDziembowski et alJl993l:lDziembowski & Pamvatnv kh| 



1993). While SPB stars pulsate with high-order g-modes 
with periods going from 1 to 3 days (frequency vA-YL /iHz),| 
low-order g and p modes as well as g-p mixed modes with 
frequencies between 30 and 100 /iHz are found excited in 
Cep stars. Miglio et al. (2007, these proceedings) show 
that, even for slow rotation, the changes of molecular gradi- 
ent induced by the turbulent mixing acting near the core lead 
to significant variations in the properties of high-order g- 
modes spectrum. Although the approximation used in Miglio| 
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et al (2007) to relate the sharpness of Nbv with the prop- 
erties of the g-mode spectra is no longer valid for low-order 
modes, the analytical description of gravity modes presented^ 
there is still able to qualitatively describe the properties of 
low-order g-mode spectra (see their Fig. 4). 

As the star evolves, the combined action of nuclear re- 
actions and convective mixing leads to a chemical composi- 
tion gradient at the boundary of the convective core, an in- 
crease of the Brunt-Vaisala frequency and therefore of the 
frequencies of gravity modes. The latter interact with pres- 
sure modes of similar frequency and affect the properties 
of non-radial oscillations by the so-called avoided crossing 
phenomenon. The modes undergoing an avoided crossing 
(mixed modes) are therefore sensitive probes of the core 
structure of the star. 

Here we study the effects of extra-mixing on the j3 Cep 
spectra by comparing the properties of low-order g and p 



for massive stars, in explaining different observational facts 
indic ating chemical mixing in the stellar radiative regions 
(e.g. iMeynet & Mae der 2005, and references therein). It is 
beyond the scope of this paper to study of the effects of rota- 
tion on the oscillation frequency spectrum. Detailed studies 
of the first, second and third order corrections of frequen- 
cies, as well as t he mode coupling d u e to rotation have been 

addressed in e.g. | Suarez et all (l2006h. Daszvriska-Daszkiewicz et aL^ 
d2002l)JSoufi et al.ldl998l) and lDziembowski & Goodd(ll992h .| 
We intend here to show that if the mechanism res ponsible 
for ch emical mixing is similar to that described in jMaederl 



2003, and references therein) additional effects will appears 



in the oscillation spectrum. In fact, Geneva models also ac- 
count for the transport of angular momentum and predict 
the rotational profile in the interior of the star as a func- 
tion of time . On the bases of the high anisotropic turbulence 
dZahnll 1992b the profile of rotational velocity generated can 



modes in models computed with overshooting and with chem-|be assumed to be shellular (f2 = Q(r)). In that case, the 
ical turbulent diffusion. Since the parameters aov and Dt adiabatic oscillation frequencies corrected by the first order 
were chosen to lead to similar evolutionary tracks, such com-| effects are given by: 
parison allows us to remove the differences in the frequen- 



cies due to a different stellar radius. In Fig. |4]we plot the 
oscillation frequencies for 10 M Q models along the main- 
sequence phase. As expected, the differences between fre- 
quencies of overshooting models and turbulent mixing ones 
is very small for pressure modes, while significant differ- 
ences appear for mixed modes. These differences increase 
with the radial order and with the angular degree of the 
modes. In the figure we have also marked with thicker sym- 
bols the frequenc ies of modes that th e non-adiabatic oscilla- 
tion code MAD jDupret et al. 2003 ) predicts to be excited. 
Here we show only the case with L>t = 7 x 10 4 cm 2 s _1 but 
computations with lower or higher efficiency of the turbu- 
lent mixing show that the differences between overshooting 
and turbulent-mixing models increase with the value of D^. 

What is most relevant from an asteroseismic point of 
view, is the change in the distance between consecutive fre- 
quencies of g-modes or modes of mixed p-g character (Az/). 
In Fig.[5]we plot these differences for the I = 2 modes of the 
10 M Q models with an effective temperature T o h ~22550 K 
(Xq — 0.3). The dots indicate differences computed be- 
tween pairs of excited modes. Therefore the difference of 
Av that we can expect between models with sharp V M (over-| 
shooting models for instance) and models with a chemi- 
cal composition gradient smoothed by the effect of, for in- 
stance, a slow rotation (V ro t ~ 50 km/s) is of the order of 
0.4 c/d (~ 5 /iHz), much larger than the precision of present 
and forthcoming observations. 

3.3 Effects of rotation 

Up to now we have considered the effect of changing V p 
by chemical mixing on the oscillation frequencies, regard- 
less of the physical processes at the origin of that chem- 
ical transport. In the last ten years, the stellar models in- 
cludin g the rotation ally induced mixing based on the the- 
ory by IZahnld 1992b have been veiy successful, particularly 



Okim = <J k m + m K ke fl(r)dr 



(1) 



where K/u are the rotation kernels based on the eigenfunc- 
tions of non-rotating model ( Lynden-Bell & Ostrikerl[l967l) 
with eigenfrequencies (7k/. 

The rotational profile for three different evolutionary stages| 
of the Geneva 10 M Q models with an initial rotational ve- 
locity Vi=50 km s _1 are shown in Fig. [6] The gradient of 
angular velocity between the convective core and the exter- 
nal layers increases with time. The ratio between angular 
velocity in the convective core and at the surface goes from 
1.3 at X c = 0.5, 1.6 at X c = 0.3, and 2.5 at X c = 0.1. 

We have estimated the first order rotational splitting for 
a uniform and a differential rotation. In the former we as- 
sume an angular velocity equal to the surface value at that 
evolutionary stage, and in the latter we use the rotational 
profile in Fig. [6] The effect of rotation d oes not significantl y 
depend on the degree i (for I < 3) (see lSuarez et alj |2006). 
In Fig. [7] we plot the rotational splitting for £ = 2 and 
\m\ = 1 modes with solid and differential rotation. Since 
the treatment of rotation evolution predicts a core rotating 
more rapidly than the surface layers, the g-modes and g- 
p mixed modes in the differential rotational case show a 
larger splitting than p-modes or modes (such as the one with 
fc=-l) whose eigenfunction has the properties of an acoustic 
one with high amplitude in the external region. On the other 
hand, the mode with k=2 is mainly concentrated in the cen- 
tral stellar regions. The rotational splitting for these modes 
ranges from 0.5 to 2 /iHz. 

There are of course also second order effects that affect 
the oscillation frequencies of rotating stars. However, for 
the low-order modes expected in (3 Cep stars, the term Vt/a 
in 50 km/s models is of the order of 2 x 10 -2 . 

For the SPB high-order modes, the Geneva 6 M Q model 
withXc = 0.3 and an initial rotation velocity Vi = 25kms~ 1 | 
provides a stellar structure with a core that rotates 1.8 times 
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Fig. 4 Frequencies of pulsation modes with angular degree t = — 3 as a function of logT c ff for main-sequence models 
of a 10 M Q star. Gray dots correspond to the frequencies of models computed with an overshooting parameter aov = 0.1, 
whereas black circles are the frequencies of models computed with a turbulent diffusion coefficient £>t = 7 x 10 4 cm s~ 2 . 
Excited modes are represented by thicker symbols. The vertical lines indicate the effective temperature of models with an 
hydrogen mass fraction at the center of the order of 0.5, 0.3 and 0. 1 (left to right). 



faster than the surface. However, given that in this model 
the rotationally induced chemical mixing is able to remove 
the sharpness in the Brunt-Vaisala frequency profile, no 
high-order g mode is particularly trapped in a confined re- 
gion near the core (see Miglio et al.2007, this proceedings) 
and therefore all the modes expected to be excited in SPB 
pulsators show similar rotational splitting. The splitting dif- 
ference between uniform and differential rotation is of the 
order of 30% and corresponds to the different angular ve- 
locity between the surface and the region where the SPB 
modes have their maximum amplitude (r/R ~ 0.3). So, the 
rotational splitting for a typical SPB is Av Iot £ 1 /iHz in 
the oscillation frequency domain going from 4 to 12 /j,Hz. 
Thus, even in the case of differential rotation we do not ex- 
pect, in the first order correction, a significantly different 
rotational splittings for different modes. However, in SPB's 
frequency domain the ratio Q/a > 0.1 and therefore second 
and third order effects of rotation may be important leading 
to asymmetric splitting of frequency modes. 



4 Conclusions 

By using a parametric modeling of turbulent mixing in the 
central region of B-type stars we have shown that the effect 
of varying the chemical gradient on the oscillation frequen- 
cies are very significant (5 /iHz for instance, for moderate 
mixing Dj- = 7 x 10 4 cm 2 s _1 ). This effect on the frequen- 
cies is in fact larger or of the same order than the rotational 
splitting for modes with a moderate rotational velocity typ- 
ical of (3 Cep (50 km s _1 ). We recall that the precision ex- 
pecte d from forthcoming observations with COROT satel- 
lite feaglin & The COROT Team! 1 19981) is of the order of 
0.1 /zH. Therefore, in asteroseismic modeling, in addition 
to the rotational splitting, the shifts of frequencies produced 
by a possible turbulent mixing should be taken into consid- 
eration. 

For high-order g-modes in SPB stars, the frequency dif- 
ferences between consecutive modes due to the trapping 
in the sharp V/i r egions is of the order of 0.08-0.7 fjHz 
(Migiio_et_al]2008). This difference almost dissapears when 
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a slight turbulent mixing (Dt = 5000 cm 2 s _1 ), that could 
result from the rotationally induced mixing in a star that has 
evolved from an initial rotation velocity of 25 km s _1 , acts 
in the central stellar region. Rotation has however other sig- 
nificant effects on SPB oscillation frequencies: a symmetric 
splitting Az; ro t > 1 fiHz due to first order corrections, and 
second and higher order corrections that can be significant 
for the low frequencies in SPBs. 

The results presented in previous sections for (3 Cep and 
SPB stars are qualitatively valid also for 5 Scuti and 7 Dor 
pulsators respectively. The situation is however even more 
complicate for those cases since these less massive stars ro- 
tate typically faster than [3 Cep and SPB. 
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Fig. 2 Upper panel: Hydrogen abundance profile in the 
central regions of 10 M models with Xq ~ 0.3. The 
different lines correspond to models computed with CLES, 
with overshooting (dashed line), and turbulent diffusion co- 
efficient £>t (solid line), and with the Geneva code with an 
initial rotational velocity of 50 km s _1 (dotted line). Lower 
panel: Brunt-Vaisala frequency profile in the central re- 
gions of the 10 M Q CLES models in upper panel. 
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Fig. 3 Total chemical diffusion coefficient for a Geneva 
10 M Q model at three different evolutionary stages, and 
with an initial rotational velocity V\ = 50 km s . The hor- 
izontal dashed line gives our chosen value for the computa- 
tions with CLES. 
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Fig. 5 Frequency difference (in c/d) between modes with 
consecutive radial order (fc) and degree I — 2 for the 10 M© 
CLES models in Fig. |2] Dots represent the differences be- 
tween theoretically predicted excited modes. 
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Fig. 6 Angular velocity profile as a function of the nor- 
malized radius for a Geneva 10 M Q model at three different 
evolutionary stages, and with an initial rotational velocity 
7i = 50kms- 1 . 




Fig. 7 First order rotational splitting for I = 2 modes of 
radial order k, for a 10 M0 model at Xc — 0.3 and with 
a uniform internal rotation rate (dashed line) and with the 
corresponding internal rotation profile give in Fig. [6] (solid 
line). 



